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Abstract: Ornithine cyclodeaminase (OCD) is an NAD + -dependent deaminase that is 
found in bacterial species such as Pseudomonas putida. Importantly, it catalyzes the direct 
conversion of the amino acid L-ornithine to L-proline. Using molecular dynamics (MD) 
and a hybrid quantum mechanics/molecular mechanics (QM/MM) method in the ONIOM 
formalism, the catalytic mechanism of OCD has been examined. The rate limiting step is 
calculated to be the initial step in the overall mechanism: hydride transfer from the 
L-ornithine' s C a -H group to the NAD + cofactor with concomitant formation of a C a =NH2 + 
Schiff base with a barrier of 90.6 kJ mol -1 . Importantly, no water is observed within the 
active site during the MD simulations suitably positioned to hydrolyze the C a =NH2 + 
intermediate to form the corresponding carbonyl. Instead, the reaction proceeds via a 
non-hydrolytic mechanism involving direct nucleophilic attack of the 5-amine at the 
C a -position. This is then followed by cleavage and loss of the 01-NH2 group to give the 
A 1 -pyrroline-2-carboxylate that is subsequently reduced to L-proline. 
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1. Introduction 

L-proline (Pro) is one of the naturally occurring 20 genetically encoded amino acids and is unique 
amongst them in being the only one with a secondary a-amine. The effect of its structure and 
functionality upon the proteins in which it is found has long been studied. In addition, however, it is 
also known to have a number of key physiologically important roles. For example, it has been shown 
to be an important antioxidant needed by microorganisms, plants, and animals [1-8]. Furthermore, in 
certain pathogenic bacteria, it has been suggested that proline metabolism plays a role in enabling the 
pathogen to survive under harsh conditions [9-11]. Plants also depend on proline production for cell 
wall biosynthesis, while mediating abiotic and biotic cell stresses [12,13]. In mammals, the 
interconversion between proline and A^pyrroline-S-carboxylate (P5C), an intermediate in its 
biosynthesis and catabolism, is believed to be involved in cell apoptosis [1,14-18]. Importantly, this 
cycle provides a redox shuttle between the cytosol and mitochondria, controlling the formation of 
reactive oxygen species [14,19,20]. Consequently, it is important to understand the functions of the 
enzymes that metabolize proline. 

In general, Pro is synthesized within cells and organisms via multi-enzymatic pathways from either 
glutamate (Glu) or arginine (Arg) [9,14,21-24]. In the former, the enzymes involved require the use of 
y-glutamyl kinase and glutamate-y-semialdehyde dehydrogenase, which require the cofactors ATP and 
NADH or NADPH, respectively. In contrast, the latter pathway utilizes the metalloenzyme arginase 
and the pyridoxal-5 f -phosphate (PLP)-dependent enzyme, ornithine 5-minotransferase. Both pathways, 
however, lead to the formation of a common intermediate, glutamic-semialdehyde. This then 
undergoes a non-enzymatic cyclization via an intramolecular condensation reaction to give the imine, 
P5C [14,25]. Then, the common enzyme P5C reductase (P5CR) reduces the latter cyclic intermediate 
to form L-Pro [9,26]. In mammals, the former pathway in which glutamate is converted to proline is 
believed to be the major metabolic route. 

More recently, however, several bacteria have been shown to be able to utilize an unusual 
enzymatic route for proline biosynthesis [27,28]. More specifically, they use the NAD + -dependent 
non-metalloenzyme ornithine cyclodeaminase (OCD) to directly convert the amino acid L-ornithine 
(L-Orn), itself also an intermediate along the above "Arg-pathway", to proline [9]. Thus, in contrast to 
the two general pathways discussed above in which the actual cyclization step is done without 
enzymatic participation, the mechanism of OCD includes the cyclization of a linear intermediate to a 
cyclic product (Scheme I) [9,29,30]. Furthermore, it not only produces L-proline stereospecifically and 
without additional enzymes being involved, but does so in an irreversible fashion [9,31]. 

Scheme I. Overall reaction for conversion of L-ornithine to L-proline as catalyzed by 
ornithine cyclodeaminase (OCD). 
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It is generally believed that the overall mechanism of OCD begins with an initial hydride 
transfer from the L-Orn substrate's C a -H moiety onto the C4 center of the nicotinamide ring of the 
NAD + cofactor. This results in the formation of an iminium (C a =NH 2 + containing) intermediate 
complex [32-37]. Unfortunately, however, it is then unclear how the mechanism proceeds. Based on 
the results of experimental studies including X-ray crystallographic structures and mass spectroscopic 
data, two possible pathways have been proposed (Scheme II). While they share some common features 
such as involving formation of a Schiff base, a chemically and biochemically important reaction 
process that has been previously studied in detail both experimentally and computationally, they also 
have some very important differences [38^7]. In particular, in the hydrolytic pathway, an active 
site water nucleophilically attacks the C a center of the iminium intermediate, resulting in loss of 
ammonium with formation of the keto acid, 5-amino-2-oxopentanoate (Scheme IIA) [28,34]. This then 
undergoes a cyclization reaction in which the 5-amine nucleophilically attacks the carbonyl carbon to 
form a cyclic A 1 -pyrroline-2-carboxylate (P2C) species. This latter intermediate has previously been 
shown to be involved in other biochemical processes including, for example, lysine catabolism. It is 
further noted that this ring formation step can be considered analogous to the cyclization of glutamic 
y-semialdehyde in which an amine (a-amine) reacts with a carbonyl carbon center. This step is then 
followed by a hydride transfer from the NADH moiety onto the exposition (C2) of P2C to 
give L-proline. 

Scheme II. Proposed (A) hydrolytic and (B) non-hydrolytic pathways for the conversion of 
L-ornithine to L-proline as catalyzed by the enzyme ornithine cyclodeaminase. 
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In contrast, in the alternate proposed mechanism known as the non-hydrolytic pathway, and shown 
above in Scheme IIB [31], the 5-amino tail of the iminium intermediate is neutralized by a suitable 
active site residue. Consequently, it is now able to nucleophilically attack the C a -center of the iminium 
intermediate to form a cyclic 2-aminoproline species. This then undergoes loss of the a-amino group as 
ammonia or ammonium to give P2C. The final step is as in the hydrolytic pathway; hydride transfer 
from the NADH moiety onto the exposition (C2) of P2C to give the final product. 

As noted above, the enzyme OCD itself has only been found in a few select bacteria [27]. However, 
it has been noted that it shows close phylogenetic resemblance to some crystallin enzymes, in 
particular those found in mammals. These proteins and enzymes are primarily known for their role in 
lens and cornea structure and transparency [48,49]. However, their malfunctioning has been linked to 
a number of diseases including cataract formation and cancer. Indeed, OCD is a member of the 
|i-crystallin family of enzymes which have been found to be abundant in the eye lens of 
marsupials [50] and are believed to have a similar chemical mechanism to that of OCD [31,51]. Thus, 
a clearer elucidation of the mechanism of OCD can provide invaluable insights not only into its 
catalytic abilities and proline biosynthesis, but also into related physiologically important enzymes. 

Computational chemistry is a proven invaluable tool for the study of enzymatic mechanisms [52]. 
Thus, we have complementarily applied molecular dynamics (MD) simulations and an ONIOM 
QM/MM approach to investigate the catalytic mechanism of OCD. In particular, we have examined 
the feasibility of both the proposed hydrolytic and non-hydrolytic pathways for conversion of 
L-ornithine to L-proline. 

2. Computational Methods 

2.1. Molecular Dynamics (MD) Equilibration 

The Molecular Operating Environment (MOE) [53] program was used for model preparation and 
the molecular dynamics (MD) simulations. For the model, the structure of OCD used was taken from 
an X-ray crystallographic structure of an OCD homo-dimer; each active site was complexed with 
L-ornithine (L-Orn) and NADH (PDB ID: 1X7D) [31]. While the protein is homo-dimeric, 
we have chosen to use a single monomer for the MD simulations. This model simplification is 
reasonable since the catalytic site does not include interface residues. Thus, a single protomeric 
enzyme-substrate-cofactor complex was selected while the other was removed. In addition, all 
selenomethionine (Se-Met) residues were mutated to the native methionines. The bound NADH 
cofactor was oxidized to give the catalytically required NAD + (i.e., a hydride was removed from the C4 
center of NADH). The coordinates of missing hydrogens were added using the MOE default method. 
A 7-A spherical layer of water molecules was then added to solvate the enzyme-substrate complex. 
In order to force the system to lie within the volume of space defined by the surrounding shell of 
waters, an ellipsoidal potential wall with a scaling constant of 2 was placed around the solvated 
enzyme-substrate complex. To allow the electrostatic and van der Waals potentials to decay smoothly, 
a damping functional factor was included. The geometry of each solvated complex was then optimized 
using the CHARMM22 force field until the root mean square gradient of the total energy fell below 
0.21 kJ mol 1 A 1 . The MD simulations were performed under constrained pressure and temperature. 
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The equations of motion were coupled with the Nose-Poincare thermostat [54] and the time step for 
numerical integration was set to 2 fs. Initially, the system was heated from 150 to 300 K for a period of 
50 ps, followed by an equilibration period of 100 ps at the constant temperature of 300 K and pressure 
of 1 atm. Based on rmsd and cluster analyses, a representative structure was chosen from the 
trajectory. This structure was then optimized with the CHARMM22 force field and used to obtain a 
suitable enzyme complex for further investigation (see below). We have successfully applied this MD 
protocol in the study of other enzymatic systems [55]. 

2.2. QM/MM Computations 

The Gaussian 09 suite of programs [56] was used for all QM/MM calculations using the ONIOM 
formalism [57-65]. Optimized structures were obtained at the ONIOM (B3LYP/6-31G(d):AMBER96) 
level of theory [66-69]. Specifically, the AMBER96 molecular mechanics force field was used to 
describe the low-layer (i.e., protein environment surrounding the active center) while the high-layer 
(i.e., the active center) was described using the density functional theory-based method B3LYP 
in conjunction with the 6-31G(d) basis set. Harmonic vibrational frequencies were calculated at 
this same level of theory in order to obtain the corresponding Gibbs free energy corrections at 
SATP (AEoibbs) in addition to characterize the nature of the stationary points on the potential 
energy surface (PES) (i.e., as energy minima or first-order transition states (TSs)). Single point 
energy calculations on the above optimized structures were performed at the ONIOM 
(B3LYP/6-311 + G(2df,p):AMBER96//B3LYP/6-31G(d):AMBER96) + AE Gi bbs level of theory within 
a mechanical embedding (ME) formalism. 

The final CHARMM22 optimized MD structure (see above) was used to obtain a suitable chemical 
model for the fully bound active site for use in the QM/MM-based mechanistic studies. More 
specifically, all residues and waters up to 15 A from the L-ornithine substrate were extracted to be used 
as the enzyme-substrate-cofactor model. This distance was chosen as it has been previously shown that 
the steric and electrostatic effects arising from the protein environment surrounding the active site 
within this distance can have important effects on the mechanism and its intermediates and transition 
states [70]. 

The QM-region included the L-ornithine substrate and the active site groups with which it directly 
interacts or are experimentally known or proposed to be involved in the mechanism. Namely, the 
R-groups of Argll2, Lys69, Glu56 and Asp228 were included along with two active site waters 
observed within the MD simulations. The NAD + cofactor was also included in part, the asterisk (*) 
indicating that only its nicotinamide and ribose ring was within the QM-region. The rest was replaced 
and modelled by a methyl group as shown in Figure 1. All those residues or waters shown in the 
outer-circle in Figure lb were included within the MM-layer in their entirety. The only exceptions 
being those in which the R-group was included within the QM region. To ensure the integrity of the 
active site model during calculations, the a-carbons of each residue as well as that of the capping 
methyl in NAD + were held fixed at their final CHARMM22 optimized positions (see above). This 
computational approach has been widely used previously, and successfully applied in studies on 
related enzymatic systems [55]. 
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Figure 1. Illustration of the ONIOM (QM/MM) fully bound active site model used: (a) the 
low or molecular mechanics (MM)-layer is shown in wire while the high or quantum 
mechanics (QM)-layer is shown in tube format; (b) those moieties within the inner circle 
were included within the QM-region, while those within the outer circle represent the 
residues contained in the MM-layer of the QM/MM model. 




3. Results and Discussion 

3.1. Structure of the Active Site with the Substrate L-Ornithine and NAD + Cofactor Bound 

We began this study by obtaining and examining the optimized structure of the fully bound active 
site complex (RC). That is, the complex in which both the NAD + cofactor and L-ornithine (L-Orn) 
substrate are bound and which is shown in Figure 2a. As it can be seen, the L-Orn is bound 
electrostatically via interactions of each of its charged groups with active site residues. In particular, 
the a-COO~ forms short, strong hydrogen bonds with the protonated R- groups of both Argll2 and 
Lys69 with distances of 1.73 and 1.53 A, respectively. Meanwhile, its a-NH 3 + group is strongly 
hydrogen bonded with the nearby R-group carboxylate of the active site aspartyl (Asp228) with an 
r(a-NH-~OOC-Asp228) distance of 1.62 A. Notably, this results in a modest lengthening in the 
participating N a -H bond to 1.08 A. For the R-group, 5-NH 3 + , however, it was found that it transferred 
a proton to the R-group carboxylate of the active site glutamyl residue (Glu56) with which it hydrogen 
bonds. Although they now form a neutral N5-HOOC-Glu56 hydrogen bond, its length is still 
decidedly short with a length of only 1.68 A. 

Importantly, in this bound conformation, the distance between the substrate's C a and the NAD + 
cofactor' s C4 center to which a hydride would be transferred, i.e., r(C a ---C4(NAD + )), is 3.98 A. It is 
noted that while in the X-ray crystal structure (PDB ID: 1X7D), the enzyme was co-crystallized with 
the substrate and NADH, the above calculated distance is in reasonable agreement with the 
experimentally measured distance of 3.8 A. It is noted that the mechanistically important C4-H-C a 
distance is 3.51 A and that the C a -H bond has lengthened marginally to 1.10 A. 
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In addition, it is noted that in the MD structure of the initial reactant complex RC and shown in 
Figure 2b, three waters (Wl, W2 and W3) were observed to lie in the active site and interact with the 
substrate. More specifically, one (Wl) forms a hydrogen bonding bridge between the 5-NH 3 + group 
and a sugar hydroxyl of the NAD + cofactor. A second (W2) is simply hydrogen bonded to the 
substrates 5-NH 3 + group while a third (W3) forms a hydrogen bonding bridge between the anionic 
R- group carboxylate of Glu56 and the substrate's protonated a-NH 3 + moiety. It should be noted that 
the oxygen of this latter water, Ow3, is 3.50 A away from the C a center. These same interactions were 
also in agreement with the subsequent QM/MM optimized structure (not shown). 

Figure 2. (a) Optimized structure (see Computational Methods) of the initial fully bound 
active site complex RC with selected bond lengths shown (in Angstroms), (b) Solvated 
MD structure of the active site. 




3.2. L-Orn Oxidation and Formation of 2-Aminoproline (AP) 

The first catalytic step of OCD involves a hydride transfer from L-Orn to NAD + via TS1 (Figure 3). 
Specifically, the hydride is transferred from the C a -H moiety of the L-Orn substrate onto the C4 center 
of the nicotinamide ring of NAD + . This process occurs via TS1 with a relatively high barrier of 
90.6 kJ mol 1 with respect to RC (Figure 4). This may reflect in part the lack of any stabilizing 
interactions between the active site and the hydride entity as well as the large structural changes 
predicted. For example, as it can be seen in Figure 3, in TS1 a significant reduction in the 
C a ---C4(NAD + ) distance is seen. Specifically, it has shortened by 1.24 A to 2.74 A. As a result, the 
hydride is almost equidistant between the C a and C4 centers with C a - • FT and C4* • H~ distances of 1.48 
and 1.28 A, respectively. However, concomitant with the FT transfer, a proton from a-NH 3 + group has 
transferred to the R-group carboxylate of Asp228 (Figure 3). Indeed, in TS1, the proton has fully 
transferred to the carboxylate as indicated by an r(H-0-Asp228) distance of 0.98 A while the C a -N 
bond length has already significantly shortened by 0.16 A to 1.35 A. This suggests that the proton 
transfer occurs early in this reaction step. 
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Figure 3. Optimized structures (see Computational Methods) with selected bond lengths 
shown (in Angstroms) for the ornithine cyclodeaminase (OCD) catalytic mechanism prior 
to deamination as L-ornithine is converted to 2-aminoproline. 




Figure 4. Potential energy surface (PES) obtained (see Computational Methods) of 
the OCD catalytic mechanism prior to deamination as L-ornithine is converted 
to 2-aminoproline. 



Relative Free Energies (kJ mol 1 ) d £ 




Complete H transfer onto the NAD cofactor results in formation of the iminium ion intermediate 
(IC1) which lies lower in energy than RC by -15.3 kJ mol" 1 (Figure 4). In IC1, the C a -C 4 (NADH) 
distance has lengthened significantly to 4.51 A and is in fact now even greater than that seen in RC 
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(Figure 3). Furthermore, the C a -N bond has shortened even further to just 1.29 A indicating the 
formation of a Schiff base. In addition, the hydrogen bond interaction between the a-NH 2 and now 
neutral Asp228-COOH group has been weakened with an r(Asp228-COOH- • -N a ) distance of 4.14 A. 

As detailed in the introduction section, it has been proposed that subsequent to hydride transfer and 
formation of the Schiff base intermediate, cyclization of L-Orn may occur via one of two possible 
pathways. Namely, it may either involve a first reaction with an active water to give the keto acid, 
5-amino-2-oxopentanoate (A20), or direct nucleophilic attack of the 5-NH2 nitrogen at the C a center. 
Indeed, the former A20 intermediate is similar to the glutamate y-semialdehyde, common in both Arg- 
and Glu-pathways for proline biosynthesis. In particular, the ketone functionality is exposed to 
nucleophilic attack by either the a-amine in the non-enzymatic reaction or by 5-amine in OCD. As 
noted previously for RC, the nearest active site water that may nucleophilically attack the C a center 
lies approximately 3.5 A from the C a center. We did attempt to elucidate a possible reaction pathway 
for reaction of IC1 with a water molecule but were unable to do so. Analysis of the structure of IC1, 
however, suggests that in agreement with Goodman et al. [31], the steric presence of Arg45 packing 
against the substrate in IC1 prevents the water from approaching much closer. Indeed, the nearest 
distance between Arg45 and the C a /N a centers is only 3.40 A. As a result, the hydrolytic pathway 
would seem unlikely to be physically feasible and hence is not being discussed further herein. Thus, 
the remainder of the discussion solely concerns the non-hydrolytic pathway unless otherwise 
mentioned. However, it is interesting to note that an iminium ion simulates a Lewis acid-activated 
carbonyl functional group [71]. In such Lewis acid activation, a lowering of the LUMO energy of the 
C=0 7r-system occurs promoting nucleophilic attack. Thus, it appears that OCD has optimized proline 
synthesis by maintaining the existence of a far more reactive intermediate such as an iminium ion 
rather than a ketone derivative. 

With the iminium intermediate IC1 formed, it is now susceptible to nucleophilic attack by N5 of 
L-Orn. However, presumably prior to attack, the "tail" is required to curl within the active site. It is 
noted that the S-NH 2 group in IC1, unlike say the a-COO~ group, forms just one hydrogen bond 
interaction with an active site residue and that being with the now neutral R- group of Glu56. This 
would be expected to facilitate the required cleavage of this interaction and subsequent curling of the 
"tail" of IC1. Unfortunately, at the present level of theory we were unable to exactly optimize a TS for 
this process. Instead, a number of detailed systematic scans of the PES for such a curling was 
performed in order to determine an upper limit for the energy required for this step. In particular, it 
was found that a decrease in the ZC a -Cp-C y -C5 by 55.4° to 128.5° gave a structure (TS2) that 
corresponded to an energy maximum of 55.2 kJ mol 1 with respect to RC, which is 70.5 kJ mol 1 with 
respect to IC1. It is noted that in TS2, the Ng-'-Ca distance has decreased from 4.51 A in IC1 to 
3.92 A. Thus, clearly the desired C a -Ns bond has not yet formed. Interestingly, it should be noted that 
all attempts to optimize energy minima complexes with ZC a -Cp-C y -C5 less than 128.5° (without 
constraining the C a -Cp-C y -C5 dihedral angle) led directly to the formation of the very low energy 
2-aminoproline intermediate complex IC2 lying 84.2 kJ mol 1 lower in energy than RC. Importantly, 
in IC2, a C a -Ns bond has been formed with a length of 1.45 A, typical for a single C-N bond. This 
suggests that in OCD, the active site structure allows for a reasonably low energy curling of the "tail" 
of the Schiff base intermediate and that once suitably positioned, the N5 center appears able to 
nucleophilically attack at the C a center readily. Thus, as stated above, given that OCD has protected 
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the iminium intermediate from reacting with water, it has allowed for an apparently barrierless C-N 
bond formation without the need for activation by a Lewis acid. This is far different than what has 
been seen in various carbonyl analogues where significant barriers for initial C-N bond formation 
exist [38-47]. 

3.3. Deamination of 2-Aminoproline (AP) 

In order to lose the a-amine of IC2 (i.e., deamination), it must first be protonated. This is analogous 
to the required protonation of the hydroxide to form water prior to Schiff base formation in the 
aldehyde analogues [38,40,41]. One possibility is that an intramolecular proton transfer may occur 
from the ring -NH 2 + - group directly onto the leaving amine. However, this would necessarily involve 
a four-membered ring transition structure which is well-known to be higher in energy due to the 
inherent strain involved [72]. Alternatively, an active site residue may be able to act as a proton donor. 
Within the active site of OCD, however, the nearest acidic residue is the neutralized R-group 
carboxylic acid of Asp228. Within the present computational model in the formation of IC2 
(i.e., during the preceding cyclization step), it was found that the Asp228-COOH group was no longer 
in a position to transfer a proton to the leaving a-amine (Figure 3). Thus, the Asp228-COOH-a-NH2 
distance must first shorten allowing, for the formation of a hydrogen bond to the leaving a-amine. 
However, for this to occur IC2 must undergo a conformational or hydrogen bonding network change. 
In particular, the proton on Glu56 must rotate allowing for a rearrangement of the active site waters 
allowing for the Asp228-COOH- • a-NH2 distance to be reduced. 

The process occurs via TS3 with a barrier of 50.6 kJ mol 1 with respect to IC2 to give the alternate 
complex IC2 f (Figure 4). The latter complex, while lying 25.4 kJ mol 1 higher in energy than IC2, 
still lies markedly lower in energy than RC by 58.8 kJ mol 1 . As expected, this rotation of the proton 
on Glu56-COOH group caused Asp228-COOH to hydrogen bond with the a-amino group of the 
substrate. Furthermore, the hydrogen bonding network between the Glu56-COOH group and 
Asp228-COOH involving the two waters also changed. Specifically, the Glu56-COOH moiety now 
acts as a hydrogen bond donor via the waters to the carboxylic group of Asp228. 

However, within the present computational model we were unable to locate a transition structure for 
proton transfer from Asp228-COOH onto the leaving amine. This is likely due to the fact that this 
would result in several charged groups (Asp228-COO~, a-NH 3 + and -NH 2 + -), all being in close 
proximity and thus highly sensitive to the degree of inclusion of the electrostatics of the protein 
environment. While it is plausible that under experimental conditions Asp228-COOH may transfer its 
proton onto the a-NH 2 group, previous computational investigations on Schiff base formation reactions 
have shown that such proton transfers generally occur via intramolecular proton transfer from the 
bridging amine; furthermore such proton transfers are facilitated by a water or other suitable group to 
prevent the formation of a four-membered TS [38,40,41]. 

Thus, we examined the possibility of a water molecule aiding proton transfer from either 
Asp228-COOH or the bridging -NH2 - to the a-NH2 moiety. Indeed, once the cyclic pyrrolidine 
species is formed in IC2 T it is noted that there is now solvent- accessible space in the vicinity of the 
Asp228-COOH group. The resulting IC2 f complex in which a water was added is hereafter denoted as 
IC2 ? -"H20 (Figure 5). The resulting PES obtained for the subsequent proton transfers is shown in 
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Figure 6a. Importantly, in IC2 ? -"H20, a water simultaneously acts as a hydrogen bond acceptor via its 
oxygen with the Asp228-COOH group and pyrrolidine's -NH 2 + - moiety, while also acting as a 
hydrogen bond donor to the a-NH 2 group of the pyrrolidine and the 2 f -OH group of the NAD + 
cofactor's ribose. Indeed, in the optimized structure of IC2 ? *"H20 r(Asp228-COOH-OH 2 ) and 
r(aN-H 2 0) were found to be 1.74 and 1.77 A, respectively, while for the ring amine r(NH 2 + -OH 2 ) 
is slightly longer at 1.95 A (Figure 5). 

Figure 5. Optimized structures (see Computational Methods) with selected bond lengths 
(in Angstroms) of the water-assisted intermediate complex (IC2 ? *"H 2 0), TS4, cyclic 
intermediate IC3, and Schiff base intermediate IC4. 




Interestingly, rather than a proton transfer from Aps228-COOH, it was found that the 
2-aminoproline itself underwent an intramolecular proton transfer from its ring -NH 2 + - moiety onto 
the a-NH 2 group via the bridging water molecule. This process occurred via the six-membered ring 
transition structure, TS4, at a very low cost of only 26.0 kJ mol -1 with respect to IC2 ? -"H 2 0 
(Figure 6a). The resulting intermediate IC3 formed lies only 13.3 kJ mol 1 higher in energy that 
IC2 ? -"H20 (Figure 6a). Importantly, the protonation of the a-amine results in a significant 
lengthening in the C a -N a bond in IC3 by 0.09 A to 1.54 A and a shortening of the C a -Ns bond 
(0.10 A) within the pyrrolidine. Such bond changes are likely to aid in deamination and 
Schiff base formation. 

In previous computational studies on Schiff base formation involving a ketone/aldehyde, the loss of 
the leaving group -OH 2 has been shown to occur in two steps [38]. First the adjacent bridging amine 
undergoes an inversion thus allowing for favorable overlap of its lone pair and the anti-bonding orbital 
of C-OH 2 bond to be cleaved. The second step is then cleavage of the C-OH 2 bond itself. For the 
2-amino-2-carboxy-pyrrolidine, no stable intermediate corresponding to inversion of the ring 
-NH- was obtained nor was a concerted TS involving cleavage of the C a -N bond with inversion of the 
bridging amine. However, it is noted that, while in previous studies, the barrier for loss of the water for 
the analogous carbonyl systems was found to be generally rate-limiting [38,39,41,42] cleavage of the 
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C a -NH 3 bond in OCD was found to be exothermic. In particular, with IC4 lying 34.9 kJ mol 1 lower 
in energy than IC3 or -21.6 kJ mol 1 relative to IC2 ? -"H 2 0 (Figure 6a). Thus, there is clearly a 
driving force for formation of this Schiff base intermediate. C-N bond cleavage with concomitant 
Schiff base formation results in the A 1 -pyrroline-2-carboxylate (P2C) containing complex IC4. It 
should also be noted that in IC4, the Asp228-COOH proton has also transferred onto the leaving NH 3 
to give an ammonium ion, NH 4 + via H 2 0. The latter is then free to leave the active site. 



Figure 6. PESs obtained (see Computational Methods) for (a) water-assisted deamination 
of 2-aminoproline to give A 1 -pyrroline-2-carboxylate (P2C) and, (b) reduction of P2C to 
give L-proline. 




3.4. Hydride Transfer from NADH onto P2C to Give the Final Product L-Proline 

For the final stage, we considered formation of the L-proline product from P2C after loss of the 
cleaved NH/ from the active site. Only minor structural changes were noticed in the resulting 
complex IC4 ? -NH4 + . In principle, L-proline can be formed via hydride transfer from the NADH 
moiety on the C2 (what will become C a ) of P2C. However, in IC4 T -NH 4 + , the key nadhCzt ■ C2 is quite 
long at 4.00 A. Furthermore, the P2C itself is not ideally positioned for the transfer. However, 
it is able to undergo an intramolecular rotation about its C2-COO" bond, i.e., a change in its 
ZOi-Ci-C a -N a dihedral angle, where the oxygen involved is the one hydrogen bonded to the Lys69. 
This process occurs via TS6 at a markedly low cost of only 14.4 kJ moF 1 with respect to IC4 ? -NH4 + 
suggesting that the rotation is likely reasonably unhindered (Figure 6b). It is noted that in TS6 the 
ZOi-Ci-C a -N a has decreased by 43.0° to 93.8°. This rotation is also exergonic with the resulting 
alternate conformer complex IC4 n being lower in energy than IC4 ? -NH4 + by 19.1 kJ mol 1 . 
Importantly, as a result of this reorientation of P2C within the active site, its C2 center is now more 
suitably positioned for hydride transfer and the nadhC4*--C2 distance has decreased significantly to 
3.14 A (Figure 7). 

Reduction of P2C via a hydride transfer from NA dhC4-H onto its C2 center can then occur via TS7 
with a barrier of 65.2 kJ mol -1 with respect to IC4" (Figure 6b). This is notably lower than the cost of 
the initial hydride transfer from the L-Orn substrate to NAD + . This final step is also exergonic with the 
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final product complex (PC) in which the L-proline is bound within the active site being a further 
8.0 kJ mol 1 lower in energy than IC4 M , and with a relative free energy of -27.1 kJ mol 1 with regards 
to IC4 ? -NH 4 + (Figure 6b). 

Figure 7. Optimized structures (see Computational Methods) with selected bond lengths 
shown (in Angstroms) of the A 1 -pyrroline-2-carboxylate intermediate IC5, TS7, and 
L-proline product (PC). 

ic4'-nh/ TS6 




4. Conclusions 

Using a combination of MD and QM/MM methods the catalytic mechanism of OCD has been 
investigated. From the results it was found that the initial hydride transfer from the C a -H group of the 
L-ornithine substrate to the C4 center of the NAD + cofactor with concomitant formation of a Schiff 
base, is the rate limiting step. In particular, this process occurred with a relative free energy barrier of 
90.6 kJ mol 1 . Experimentally, this hydride transfer step might be examined by deuterating 
H-(C a ) of L-Orn and conducting a kinetic isotope effect study. 

For the remaining steps of the mechanism, while two pathways have been proposed, it appears that 
the enzyme most likely operates via a non-hydrolytic pathway. In particular, the MD and QM/MM 
results suggest that water is sterically hindered from attacking C a after the initial Schiff base formation. 
Given that iminium ions are generally more reactive it seems that the active site of OCD evolved to 
guarantee such an intermediate exists by preventing its reaction with water. Indeed, following a 
conformational change of the substrate within the active site, a barrierless C-N bond formation 
occurred. This is considerably different than seen in previous investigations of various carbonyl 
analogues where significant barriers to C-N bond formation exists. While a transition state for 
deamination could not be found, the overall process was found to be thermodynamically favorable. 
Importantly, with deamination, a second Schiff base was formed. Like the initial C-N bond formation, 
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this Schiff base would likely become more reactive. Indeed, the final step in the reaction was a FT 
transfer with a low barrier of 65.2 kJ mol 1 . The resulting product, L-proline, being thermodynamically 
more favored than the preceding intermediate. 

Thus, while the typical biosynthetic pathway for L-proline from arginine requires two enzymes and 
non-enzymatic cyclization of a glutamate y-semialdehyde intermediate, this cyclic intermediate being 
common to the pathway involving glutamate, is later reduced to form L-Pro. In contrast, ornithine 
cyclodeaminase appears to first exploit the formation of a highly reactive C a =NH2 + -containing 
iminium ion in order to enable cyclization and ultimately formation of a second iminium ion, 
A 1 -pyrroline-2-carboxylate, formed after loss of the a-NH 3 group. The latter ion is readily reduced via 
a hydride transfer from the NADH cofactor onto its "C a -center", resulting in the formation 
of L-proline. 
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